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Abstract 

The shadow of a rotating black hole with nonvan- 
ishing gravitomagnetic charge has been studied. It was 
shown that in addition to the angular momentum of 
black hole the gravitomagnetic charge term deforms the 
shape of the black hole shadow. From the numerical 
results we have obtained that for a given value of the 
rotation parameter, the presence of a gravitomagnetic 
charge enlarges the shadow and reduces its deformation 
with respect to the spacetime without gravitomagnetic 
charge. Finally we have studied the capture cross sec- 
tion for massive particles by black hole with the nonva- 
nishing gravitomagnetic charge. 

Keywords Photon motion Shadow of Black hole NUT 
spacetime 



1 Introduction 

At present there is no any observational evidence for 
the existence of gravitomagnetic monopole, i.e. so- 
called NUT ([Newman et al.lll963l ) parameter. There- 
fore study the motion of massless and the massive 
test particles and particle acceleration mechanisms in 
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NUT spacetime may provide new tool for studying 
new important general relativistic effects which are 
associated with nondiagonal components of the met- 
ric tensor and have n o Newtonian analogu e s (See , 
e^ 



Nour-Zonoz (2004) 



Morozova and Ahmedovl (|2009t ) where 



Kagramanova et al 



solutions for 

electromagnetic waves and interferometry in spacetime 
with NUT parameter have been studied.). Kerr-Taub- 
NUT spaceti me with Maxwell a nd dilation fields is in- 
vesti gated by Aliev et al. (|2008h ■ In our preceding pa - 



pers (jMorozova et al.ll2008t lAbduiabbarov et al.l l2008) 



we have studied the plasma magnetospherc around a 
rotating, magnetized neutron star and charged particle 
motion around compact objects immersed in external 
magnetic field in the presence of the NUT parame- 
ter. The Penrose process in the spacetime of rotating 
black hole with nonva nishing gravitomagnetic ch arge 
has been considered by lAbduiabbarov et al.l (|201ll ). 

Analytic treatment of geodesies of electrically and 
magnetically charged test particles in the Reissner- 

Nordstro m and Taub-NUT space-times hav e been con- 

sidere d by Grunau and Kagramanova ( 2011 ) ; Kagramanova et aL^ 
( 2010tk Hackmann and Lammerzahl ( 2012 ). Kubiznak 
and coauthors have considered the parallel-propagated 
frame along null geodesies in higher-dim ensional black 



hole spacetimes ( Kubiznak et al. 20091) . The parallel 



transport equations in the higher-dimen sional Kerr 
NUT-( A)dS spacetimes hav e been solved by Cormel et alj 
( 20081 ) : lFroiov and Krtousl (|201lh : lKrtous et al.l (|2008l ) 
Magnetized black hole on the Taub-NUT instan- 
ton has been considered by the authors of the pa- 
per ( Nedkova and Yazadiievll2012 ). Some geometrical 
properties of Taub-NUT metric, such as asymptotically 
fla tness and varia tional principle have been analyzed 
bv lVirmanil |201ll ). 

In this paper we study photon orbits around the ro- 
tating black hole with nonvanishing gravitomagnetic 
charge as well as black hole shadow in the vicinity 



2 



of Kcrr-Taub-NUT spacetime. The paper organized 
as follows: in Sec. 2, we review the basic aspects of 
the geometry and the geodesies of the Kerr-Taub-NUT 
black hole. In Sec. 3, we obtain the shadows of black 
holes for the different values of the angular momentum 
and gravitomagnetic charge of black hole. In Sect. 4 
we study the capture cross section for massive parti- 
cles by black hole with nonvanishing gravitomagnetic 
charge. Finally, in Sec. 5, we discuss the results found. 
Throughout the paper, we use a space-like signature 
(— , +, +, +) and a system of units in which G = 1 = c 
(However, for those expressions with an astrophysical 
application we have written the speed of light explic- 
itly.). Greek indices are taken to run from to 3 and 
Latin indices from 1 to 3; covariant derivatives are de- 
noted with a semi-colon and partial derivatives with a 
comma. 



2 Photon motion around Kerr-Taub NUT 
black holes 

Here we will study massless particles motion in the 
vicinity of a black hole of mass M in the presence 



time metric (Newman et al.ll963t 


Morozova et al.ll2008 


Abdujabbarov et al. 2008): 



ds 2 = - i (A - a 2 sin 2 9) dt 2 + Z + dd 2 



[(£ + ax ) 2 sin 2 0- X 2 A] 



+— (Ax - a(E + a\) sin 2 6)dtd<j>, 
where notations A, E, and x are denned as 



(i) 



A = r 2 +a 2 - 2Mr~l 2 , 



X 



! sin 2 9 — 21 cos 9, 



here a is the specific angular momentum per total mass 
of black hole M(a = J/M) and I is the gravitomagnetic 
charge. The event horizon is determined by largest root 
of the equation A = 0, given by 



M + ^M 2 -a? + l 2 . 



(2) 



Consider a black hole placed between a source of 
light and an observer. Then the light reaches the ob- 
server after being deflected by the black hole gravita- 
tional field. But some part of the deflected light with 



small impact parameters emitted by the source falling 
into the black hole. This phenomena result a dark fig- 
ure in the map of the space called the shadow. The 
boundary of this shadow defines t he sh ape of a black 
hole (see e.g. lAmarilla and Eiroal (|2012|)). The study 
of the geodesic structure around black hole is very im- 
portant to obtain the apparent shape. The Hamilton- 
Jacobi equation determines the geodesies for a given 
space-time geometry: 



dS_ 
~dX 



1 



, dS dS 
dx 11 dx v ' 



(3) 



where A is an affine parameter along the geodesies, 
g^ v are the components of the metric tensor and 
S is the Jacobi action. If the problem is separa- 
ble (the seperable problem of Hamilton- Jacobi equa- 
tio n in Kerr-Taub-NUT spacetim e has been studied 
by iDadhich and Turakulovl ( 2002h ). the Jacobi action 
S can be written in the form 



3 = ~m 2 A -£t + £<f> + S r {r) + S g {9) , 



(4) 



where m is the mass of a test particle. The second term 
in the right hand side is related to the conservation of 
energy £, while the third term is related to the conser- 
vation of the angular momentum C in the direction of 
the axis of symmetry. In the case of null geodesies, we 
have that m — 0, and from the Hamilton- Jacobi equa- 
tion, the following equations of motion are obtained: 



dX 



E dA 



+ a 2 + l 2 



[('< 



l 2 )£ - aC] 



sin 9 
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dr 
dX 

de 

dX 



sin 

n, 



+ l 2 )£ 



aC 



(5) 

(6) 
(7) 
(8) 



where the functions lZ(r) and <d(9) are defined by 



K = 
6 = 



[(» 



K, + cos 2 9 



+4l£ cost 



l 2 )£- 
(a 2 - 

£a — 



aC\' 
Al 2 

sin 2 t 
C 

sin 2 6 



A [AC + (£ 
C 2 



)£*- 



sm 



a£) 2 \<)) 



(10) 



and /C is a constant of separation. Defining the effective 
potential for massless particles as (dr/dX) 2 = V c g one 
may study the radial motion of photons in the presence 
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of gravitomagnetic charge. In the Fig. [T] the radial 
dependence of the effective potential of radial photon 
motion is shown. From the figure it is seen that with 
the incease of the gravitomagnetic charge the shape of 
the effective potential is going to shift to the observer 
at infinity. This corresponds to increasing the event 
horizon of the Kerr-Taub-NUT black hole. Moreover, 
one may conclude from the Fig. [T]that with the increase 
of the gravitomagnetic charge the circular photon orbits 
become unstable. 



0.5 




Fig. 1 The radial dependence of the effective potential 
of radial motion of the massless particles for the different 
values of the gravitomagnetic charge: solid line for l/M = 
0.1, dashed line for l/M = 0.5, and dot-dashed line for I = 
0.9. 

The propagation of light in the Kerr-Taub-NUT 
spacetime is totally determined by the equations ([S"])- 
(|8]l . Two impact parameters can be expressed in terms 
of the constants of motion £, C and the Carter con- 
stant JC and characterize the light rays. One can easily 
define impact parameters for general orbits around the 
black hole as £ = C/£ and rj — K,/£ 2 . The equation 
can be used to derive the orbits with constant r 
in order to obtain the boundary of the shadow of the 
black hole. These orbits satisfy the conditions fulfilled 
by the values of the impact parameters that determine 
the contour of the shadow lZ(r) = = dlZ(r)/dr, , 

= a 2 (l + r)+r 2 (r-3) + Z 2 (l-3^ 
a(l — r) 

r)(r) = a- 2 (r-l)- 2 jr 3 [4a 2 -r(r-3) 2 ]-/ 2 [4r 2 a 2 



+ (1 - 3r)(/ 2 (l - 3r) - 6r 2 + 4a 2 r + 2r 3 )] j (.12) 

Note that the corresponding values being valid in the 
Kerr geometry can be obtained in the limiting case 
when I = 0. 



3 Kerr-Taub-NUT black hole shadow 



observer 




Fig. 2 The schematic geometry of the gravitational lens. 
An observer far away from the black hole, can set up a ref- 
erence coordinate system (x, y, z) with the black hole at the 
origin. The Boyer-Lidquist coordinates coincide with this 
system only at infinity. The tangent vector to an incoming 
light ray defines a straight line, which intersects the a — ft 
plane at the point (ai,/3i). 

Using the celestial coordinate s one can easily de- 
scribe the shadow (see for example Vazquez and Estebanfr 
( 2004 ): 



lim I — rn smfo 
r ->oo \ dr 



and 



ft = hm r — 

r ->oo dr 



(13) 



(14) 



since here an observer far away from the black hole is 
considered r — > oo, is the angular coordinate of the 
observer, i.e. the inclination angle between the rota- 
tion axis of the black hole and the line of sight of the 
observer. The geometry of the new introduced coordi- 
nates is shown in Fig. [5] The coordinates a and ft are 
the apparent perpendicular distances of the image as 
seen from the axis of symmetry and from its projection 
on the equatorial plane, respectively. 

Calculating d<j)/dr and d9/dr using the spacetime 
metric (JT]) and taking the limit of a far away observer 
one can obtain celestial coordinates as 



a = — £csc( 



(15) 



•1 



and 

p = ± 



and 



77 + (a 2 



4/ 



cos 2 6> -£ 2 cot 2 6»o 



sm fg 



+4Z cos 6> (a 



sin 2 6*o ' 



1/2 



(16) 



where expressions (j6|), ([7]), and ([8]) were used to calcu- 
late dO/dr and dip /dr. These equations have implicitly 
the same form as for the Kerr metric, with the new ra- 
dial functions £ and r\ given by Eqs. (jlip and (|T2|) (a 
detailed calculation of the values of £ and -q, and the 
expressions of the celestial coordinates a and j3 as a 
function of the co nstants of motion for the Ke rr geom- 
etry, are given in ( Vazquez and Esteban 20041 )). 

In the case of rotating black hole, one may introduce 
two observables which approximately characterize the 
apparent shape. First one should approximate the ap- 
parent shape by a circle passing through three points 
which are located at the top position, bottom, and the 
most right end of the shadow as shown in Fig. [3] the 
radius R s of the shadow is defined by the radius of this 
circle. One can also define the distortion parameter S s 
of the black hole shadow by S s = D cs /R s . Two vari- 
ables (R s and 8 S ) can be interpreted as observab le in 
astronomical observation ( Hioki and Maedall2009t) . 




(18) 



For the visualization of the shape of the black hole 
shadow one needs to plot f3 vs a. In Fig. HI we show the 
contour of the shadows of black holes. From the plots 
one can see that the presence of the gravitomagnctic 
charge will increase the effective size of the shadow. 
We plot the shapes of the shadow of black hole with the 
gravitomagnetic charge for the different values of black 
hole angular momentum a: a/M — 0.5, a/M = 0.7, 
a/M = 0.8, and a/M = 0.99. One can compare the 
effect of the NUT parameter and the black hole angular 
momentum on modification of the shape of the shadow 
of black hole. It appears they have opposite effects 
on black hole shadow size. The gravitomagnetic charge 
increases the size of the shadow shape while black hole's 
angular momentum decreases its size. 

The observable R s can be calculated from the equa- 
tion 



(a t - a,.) 2 + fl 2 t 



and the observable S s is given by 



R s ' 

where (a p ,0) and (a p ,0) are the points where the ref- 
erence circle and the contour of the shadow cut the 
horizontal axis at the opposite side of (cv r ,0), respec- 
tively. In Fig. [5J the observables R s and S s are shown 
as functions of the gravitomagnetic charge /. From the 
dependence of R s from the NUT parameter one can 
again see that the gravitomagnetic charge forces to in- 
crease the size of the black hole shadow. The depen- 
dence of 8 S from NUT charge shows that gravitomag- 
netic charge forces to shadow to get the shape of circle 
than ellipse. In the case of rotation, with the increase 
of black hole's angular momentum the shape of black 
hole shadow takes form of ellipse rather than circle. 



Fig. 3 The observables for the apparent shape of a rotating 
black hole are the radius R s and the distortion parameter 
Ss = D cs /R s . Here D cs is the difference between the left 
endpoints of the circle and of the shadow. 

If the observer is situated in the equatorial plane of 
the black hole, then the inclination angle is 8q — ir/2. 
The gravitational effects on the shadow, which grow 
with 8q, are larger. In this case, one has 



4 PARTICLE CAPTURE CROSS SECTIONS 
FOR BLACK HOLE WITH NUT CHARGE 

In this section we will study the pure effect of NUT pa- 
rameter assuming that the the angular momentum of 
the black ho le is equal to zero. It has been shown in 



(17) 



the paper of lAbdujabbarov et al.1 (|20081 ) that variables 
in the Hamilton- Jacobi equation for the particle motion 
around NUT black hole can be separated in the equa- 
torial plane. In the space-time metric ([T]) we assume 
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that the central object is non-rotating and particles are 
confined at the equator ial plane (a = 0, and 8 = tt/2) . 

It was first shown by IZimmerman and Shahir ( 19891 ) 
for the spherical sym metric case (NUT s pacetime) and 
later in the paper of Bini et aT] ( 20031) for the axial 
symmetric case (Kerr-Taub-NUT spacetime) that the 
orbits of the test particles are confined to a cone with 
the opening angle 8 given by cos 8 = 2£l/C. It also 
follows that in this case the equatio ns of motion on 
the cone depend on I o nly via I 2 ( Bini et ahl 2003 ; 
Abduiabbarov et al.ll2008l ). 

The main point is that the small value for the up- 
per limit for gravitomagnetic moment has been ob- 
tained by comparing theoretical results with experi- 
mental data a s (i) l/M < 10~ 24 from th e gravitational 
microlensing (jRahvar and Habibil 120041 ) . (ii) l/M < 
1.5-10 -18 fro m the interferometry experimen ts on ultra- 
cold atoms ( Morozova and Ahmedov 20091 ), (hi) and 
similar limit has been obtained fr om the experiments where 
on M ach-Zchnder interferometer (jKagramanova et al. 
120081) (here M is the total mass of central gravitating p 
object). 

Due to the smallness of the gravitomagnetic charge 
let us consider the motion in the quasi-equatorial plane 
when the motion in 8 direction changes as 8 = tt/2 + 
68(t), where 88(t) is the term of first order in I, then it 
is easy to expand the trigonometric functions as sin 8 = 
1 - S8 2 (t)/2 + 0(S8 4 (t)) and cosfl = 58(t) - 0(S8 3 (t)). 
Neglecting the small terms 0(S8 2 (t)), one can write 

£ = r 2 + I 2 , A = r 2 - 2Mr - I 2 , and \ = an( i con - C? CT = 16M 2 - 6/ 2 
scquently the equation of motion for the radial motion 
takes the following form 



r/M 

Fig. 6 The radial dependence of the effective potential of 
radial motion of the massive particles for the different values 
of the gravitomagnetic charge: solid line for l/M = O.f, 
dashed line for l/M — 0.5, and dot-dashed line for I = 0.9. 



r 



/ = 



C = 



Gravitational capture of the particle occurs for C < 
£ CT . For the L = C cr orbit spirals into a circular orbit, 
the radius of which is determined by the value of the 
multiple root of the polynomial (f2Tj) . i.e. discriminant 
of the later should vanish. Now it is easy to find the 
expression for £ cr as 



13Z 4 
16M 



(22) 



R(r) - [£ 2 -l-2U oB (r,l,C)] 



(19) 



where £ and C is the energy and angular momentum of 
the particle per unit of its mass and the quantity 



U eS {r, l,C) = - 



I 2 + Mr AC 2 



2Y? 



(20) 



can be interpreted as effective potential of the radial 
motion of the test particle at equatorial plane. The 
radial dependence of the effective potential of radial 
motion of the massive particles for the different values 
of the gravitomagnetic charge is presented in Fig. [6j 

Assuming that the uncharged particle is moving 
slowly at infinity, i.e. £ ~ 1 one can easily rewrite 
the expression (|19j) in the following form: 




p 3 + (f 



C)p 2 + C P 



IC 
~2 



(21) 



Fig. 7 The dependence of the critical angular momentum 
for capturing by central black hole from gravitomagnetic 
monopole momentum. 

In the FigJ7]the dependence of C CT from dimension- 
less NUT parameter is presented. The dependence 
shows that the presence of the gravitomagnetic charge 
decreases the capture cross section for particles by black 
hole. 
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5 Conclusion 

In this paper, we have studied the shadow of black 
hole with nonvanishing gravitomagnetic charge and an- 
alyzed how the shadow of the black hole will be dis- 
torted by the presence of the NUT parameter. From 
the numerical results we have obtained that the NUT 
parameter forces to increase the size of the black hole 
shadow. The dependence of the distortion parameter S s 
from the NUT charge shows that the gravitomagnetic 
charge forces black hole's shadow to get the shape of cir- 
cle than ellipse. We have also studied the capture cross 
section for massive particles by black hole with non- 
vanishing gravitomagnetic charge and found its strong 
dependence from the NUT parameter. 
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Fig. 4 Silhouette of the shadow cast by a Kerr-Taub-NUT black hole situated at the origin of coordinates with inclination 
angle 9 — it/2, having a black hole angular momentum a and a NUT charge I. Upper row, left: a/M — 0.5, l/M — 0.1 
(solid line), l/M = 0.5 (dashed line), and l/M = 0.9 (dashed-dotted line). Upper row, right: a/M = 0.7, l/M = 0.1 (solid 
line), l/M = 0.5 (dashed line), and l/M = 0.9 (dashed-dotted line). Lower row, left: a/M = 0.8, l/M = 0.1 (solid line), 
l/M = 0.5 (dashed line), and l/M = 0.9 (dashed-dotted line). Lower row, right: a/M = 0.99, l/M = 0.1 (solid line), 
l/M — 0.5 (dashed line), and l/M = 0.9 (dashed-dotted line). The shadow corresponds to each curve and the region inside 
it. 
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Fig. 5 Observables R a and 8 a as functions of the NUT charge, corresponding to the shadow of a black hole situated at 
the origin of coordinates with inclination angle 9 — tt/2 and spin parameters a = 0.99. 
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Abstract 

The shadow of a rotating black hole with nonvan- 
ishing gravitomagnetic charge has been studied. It was 
shown that in addition to the angular momentum of 
black hole the gravitomagnetic charge term deforms the 
shape of the black hole shadow. From the numerical 
results we have obtained that for a given value of the 
rotation parameter, the presence of a gravitomagnetic 
charge enlarges the shadow and reduces its deformation 
with respect to the spacetime without gravitomagnetic 
charge. Finally we have studied the capture cross sec- 
tion for massive particles by black hole with the nonva- 
nishing gravitomagnetic charge. 

Keywords Photon motion Shadow of Black hole NUT 
spacetime 



1 Introduction 

At present there is no any observational evidence for 
the existence of gravitomagnetic monopole, i.e. so- 
called NUT ([Newman et al.lll963l ) parameter. There- 
fore study the motion of massless and the massive 
test particles and particle acceleration mechanisms in 
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NUT spacetime may provide new tool for studying 
new important general relativistic effects which are 
associated with nondiagonal components of the met- 
ric tensor and have n o Newtonian analogu e s (See , 
e^ 



Nour-Zonoz (2004) 



Morozova and Ahmedov! (|2009t ) where 



Kagramanova et al 
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solutions for 

electromagnetic waves and interferometry in spacetime 
with NUT parameter have been studied.). Kerr-Taub- 
NUT spaceti me with Maxwell a nd dilation fields is in- 
vesti gated by Aliev et al. (|2008h ■ In our preceding pa - 



pers (jMorozova et al.ll2008t lAbduiabbarov et al.l l2008) 



we have studied the plasma magnetospherc around a 
rotating, magnetized neutron star and charged particle 
motion around compact objects immersed in external 
magnetic held in the presence of the NUT parame- 
ter. The Penrose process in the spacetime of rotating 
black hole with nonva nishing gravitomagnetic ch arge 
has been considered by lAbduiabbarov et al.l (|201ll ). 

Analytic treatment of geodesies of electrically and 
magnetically charged test particles in the Reissner- 

Nordstro m and Taub-NUT space-times hav e been con- 

sidere d by Grunau and Kagramanova ( 2011 ) ; Kagramanova et aL^ 
( 2010tk Hackmann and Lammerzahl ( 2012 ). Kubiznak 
and coauthors have considered the parallel-propagated 
frame along null geodesies in higher-dim ensional black 



hole spacetimes ( Kubiznak et al. 20091) . The parallel 



transport equations in the higher-dimen sional Kerr 
NUT-( A)dS spacetimes hav e been solved by Cormel et alj 
( 20081 ) : lFroiov and Krtousl (|201lh : lKrtous et al.l (|2008l ) 
Magnetized black hole on the Taub-NUT instan- 
ton has been considered by the authors of the pa- 
per ( Nedkova and Yazadiievll2012 ). Some geometrical 
properties of Taub-NUT metric, such as asymptotically 
fla tness and varia tional principle have been analyzed 
bv lVirmanil |201ll ). 

In this paper we study photon orbits around the ro- 
tating black hole with nonvanishing gravitomagnetic 
charge as well as black hole shadow in the vicinity 



2 



of Kcrr-Taub-NUT spacetime. The paper organized 
as follows: in Sec. 2, we review the basic aspects of 
the geometry and the geodesies of the Kerr-Taub-NUT 
black hole. In Sec. 3, we obtain the shadows of black 
holes for the different values of the angular momentum 
and gravitomagnetic charge of black hole. In Sect. 4 
we study the capture cross section for massive parti- 
cles by black hole with nonvanishing gravitomagnetic 
charge. Finally, in Sec. 5, we discuss the results found. 
Throughout the paper, we use a space-like signature 
(— , +, +, +) and a system of units in which G = 1 = c 
(However, for those expressions with an astrophysical 
application we have written the speed of light explic- 
itly.). Greek indices are taken to run from to 3 and 
Latin indices from 1 to 3; covariant derivatives are de- 
noted with a semi-colon and partial derivatives with a 
comma. 



2 Photon motion around Kerr-Taub NUT 
black holes 

Here we will study massless particles motion in the 
vicinity of a black hole of mass M in the presence 



time metric (Newman et al.ll963t 


Morozova et al.ll2008 


Abdujabbarov et al. 2008): 



ds 2 = - i (A - a 2 sin 2 9) dt 2 + Z + dd 2 



[(£ + ax ) 2 sin 2 0- X 2 A] 



+— (Ax - a(E + a\) sin 2 6)dtd<j>, 
where notations A, E, and x are denned as 



(i) 



A = r 2 +a 2 - 2Mr~l 2 , 



X 



! sin 2 9 — 21 cos 9, 



here a is the specific angular momentum per total mass 
of black hole M(a = J/M) and I is the gravitomagnetic 
charge. The event horizon is determined by largest root 
of the equation A = 0, given by 



M + ^AI 2 -a 2 + l 2 . 



(2) 



When a black hole is placed between a source of light 
and an observer, the light reaches the observer after be- 
ing deflected by the black hole gravitational field; but 
some part of the photons emitted by the source, those 



with small impact parameters, end up falling into the 
black hole, not reaching the observer, giving as result a 
dark zone in the sky called the shadow. The apparent 
shape of a black hol e is thus defined by the bou ndary of 
the shadow (see e.g. lAmarilla and Eiroal (|2012h ). In or- 
der to obtain the apparent shape, we need to study the 
geodesic structure around black hole. The Hamilton- 
Jacobi equation determines the geodesies for a given 
space-time geometry: 



dS_ 
~dX 



1 



, dS dS 
dx 11 dx v ' 



(3) 



where A is an affine parameter along the geodesies, 
g^ v are the components of the metric tensor and 
S is the Jacobi action. If the problem is separa- 
ble (the seperable problem of Hamilton- Jacobi equa- 
tio n in Kerr-Taub-NUT spacetim e has been studied 
by iDadhich and Turakulovl ( 2002h ). the Jacobi action 
S can be written in the form 



3 = ^m 2 X -£t + £<f> + S r {r) + S g {9) , 



(4) 



where m is the mass of a test particle. The second term 
in the right hand side is related to the conservation of 
energy £, while the third term is related to the conser- 
vation of the angular momentum C in the direction of 
the axis of symmetry. In the case of null geodesies, we 
have that m — 0, and from the Hamilton- Jacobi equa- 
tion, the following equations of motion are obtained: 



dX 



+ a 2 + l 2 



[('< 



l 2 )£ - aC] 



sin 9 



A 

+ 



1 



dr 
dX 

de 

dX 



sin 

n, 



+ l 2 )£ 



aC 



(5) 

(6) 
(7) 
(8) 



where the functions lZ(r) and <d(9) are defined by 



K = 
6 = 



[(» 



K, + cos 2 9 



+4l£ cost 



l 2 )£- 
(a 2 - 

£a — 



aC\' 
Al 2 

sin 2 t 
C 

sin 2 6 



A [AC + (£ 
C 2 



)£*- 



sm 



a£) 2 \<)) 



(10) 



and /C is a constant of separation. Defining the effective 
potential for massless particles as (dr/dX) 2 = V c g one 
may study the radial motion of photons in the presence 



3 



of gravitomagnetic charge. In the Fig. [T] the radial 
dependence of the effective potential of radial photon 
motion is shown. From the figure it is seen that with 
the incease of the gravitomagnetic charge the shape of 
the effective potential is going to shift to the observer 
at infinity. This corresponds to increasing the event 
horizon of the Kerr-Taub-NUT black hole. Moreover, 
one may conclude from the Fig. [T]that with the increase 
of the gravitomagnetic charge the circular photon orbits 
become unstable. 



0.5 




Fig. 1 The radial dependence of the effective potential 
of radial motion of the massless particles for the different 
values of the gravitomagnetic charge: solid line for l/M = 
0.1, dashed line for l/M = 0.5, and dot-dashed line for I = 
0.9. 
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Fig. 2 The geometry of the gravitational lens. An observer 
far away from the black hole, can set up a reference coordi- 
nate system (x, y, z) with the black hole at the origin. The 
Boyer-Lidquist coordinates coincide with this system only 
at infinity. The reference frame is chosen so that, as seen 
from infinity, the black hole is rotating around the z axis. 
In this system, the line joining the origin with the observer 
is normal to the a — /3 plane. The tangent vector to an in- 
coming light ray defines a straight line, which intersects the 
a — /3 plane at the point (<Xi, j3i). 



The equations ©-(JSJ) determine the propagation of 
light in the Kerr-Taub-NUT spacetime. The light rays 
are, in general, characterized by two impact parame- 
ters, which can be expressed in terms of the constants 
of motion £ , C and the Carter constant /C. Combin- 
ing these quantities we define as usual £ = C/£ and 
r\ — JC/S 2 , which are the impact parameters for gen- 
eral orbits around the black hole. We use Eq. to 
derive the orbits with constant r in order to obtain the 
boundary of the shadow of the black hole. These orbits 
satisfy the conditions 72, (r) = = dlZ(r)/dr, which are 
fulfilled by the values of the impact parameters that 
determine the contour of the shadow, namely, 



£(r) 



a 2 (l + r)+r 2 {r- 3) + Z 2 (l-3r) 
a(l-r) : 

>- 2 fr 3 [4a 2 -r(r-3) 2 ]- 



(11) 



a-\r 



l 2 Ur 2 a 2 



+ (1 - 3r)(7 2 (l - 3r) - 6r 2 + 4a 2 r + 2r 3 )] \ (.12) 



Note that the corresponding values being valid in the 
Kerr geometry can be obtained in the limiting case 
when I = 0. 



Adopting the celestial coordinates is very convenient 

to des cribe the shadow (see for example Vazquez and Esteb anj 
( 2004h h 



lim 



and 



lim 



-r sin 6»o 



d6 

dr 



dr 



(13) 



(14) 



since we consider an observer far away from the black 
hole ro goes to infinity, 9q is the angular coordinate of 
the observer, i.e. the inclination angle between the ro- 
tation axis of the black hole and the line of sight of the 
observer. The geometry of the new introduced coordi- 
nates is shown in Fig. [2] The coordinates a and /3 are 
the apparent perpendicular distances of the image as 
seen from the axis of symmetry and from its projection 
on the equatorial plane, respectively. 

Calculating d<f>/dr and dO/dr from the space-time 
metric given by expression (fT|) and taking the limit of 
a far away observer we obtain celestial coordinates as a 
function of the constants of motion in the form 



a = — £ esc ( 



(15) 



•1 



and 

-- 



± 



T)+ (a 2 



4/ 



, ; cos 2 6 Q - £ 2 cot 2 6> 



+4Z cos 6*o (a 



sm t>o 



sin 2 # 



1/2 



(16) 



where Eqs. ((6|), (jTj) , and (0 were used to calculate 
dB I dr and d(j)/dr. These equations have implicitly the 
same form as for the Kerr metric, with the new ra- 
dial functions £ and r\ given by Eqs. and (|T2"j) (a 
detailed calculation of the values of £ and -q, and the 
expressions of the celestial coordinates a and f3 as a 
function of the co nstants of motion for the Ke rr geom- 
etry, are given in (jVazquez and Estebanll2004l )). 

In the case of a Kerr black hole, we may introduce 
two observables which approximately characterize the 
apparent shape. First we approximate the apparent 
shape by a circle passing through three points which 
are located at the top position (A), the bottom posi- 
tion (B), and the most right end (C) of the shadow 
as shown by three red points in Fig. [3] The point C 
corresponds to the unstable retrograde circular orbit 
when seen from an observer on the equatorial plane. 
We define the radius R s of the shadow by the radius of 
this approximated circle. We also take into account the 
dent in the left-hand side of the shadow. The size of 
this dent is evaluated by D cs , which is the difference be- 
tween the left endpoints of the circle and of the shadow 
(see Fig. 3). Then we define the distortion parameter 
S s of the shadow by S s = D cs /R s . Thus we adopt these 
two variable s (R* and ) as observa ble in astronomical 
observation ([Hioki and Maedal l2009'). 

When the observer is situated in the equatorial plane 
of the black hole, the inclination angle is do = tt/2 
and the gravitational effects on the shadow, which grow 
with #o, are larger. The inclination angle corresponding 
to the Galactic supermassive black hole is also expected 
to lie close to n/2. In this interesting case, we have 
simply 



a = -£ , 



and 



(17) 



(18) 



For the visualization of the shape of the black hole 
shadow one needs to plot f3 vs a. In Fig. 2J we show the 
contour of the shadows of black holes. From the plots 
one can see that the presence of the gravitomagnetic 
charge will increase the effective size of the shadow. 
We plot the shapes of the shadow of black hole with the 
gravitomagnetic charge for the different values of black 




Fig. 3 The observables for the apparent shape of a Kerr 
black hole are the radius R 3 and the distortion parameter 
5 S — D cs /R a , approximating it by a distorted circle, where 
D C3 is the difference between the left endpoints of the circle 
and of the shadow. 



hole angular momentum a: a/M = 0.5, a/M = 0.7, 
a/M = 0.8, and a/M = 0.99. One can compare the 
effect of the NUT parameter and the black hole angular 
momentum on modification of the shape of the shadow 
of black hole. It appears they have opposite effects 
on black hole shadow size. The gravitomagnetic charge 
increases the size of the shadow shape while black hole's 
angular momentum decreases its size. 

The observable R s can be calculated from the equa- 
tion 



R, = 



(a t ~a r ) 2 + f3 2 
2\a t - a r \ 



and the observable S s is given by 



R s ' 

where (<S P ,0) and (a p ,0) are the points where the ref- 
erence circle and the contour of the shadow cut the 
horizontal axis at the opposite side of (cv r ,0), respec- 
tively. In Fig. the observables R s and S s are shown 
as functions of the gravitomagnetic charge I. From the 
dependence of R s from the NUT parameter one can 
again see that the gravitomagnetic charge forces to in- 
crease the size of the black hole shadow. The depen- 
dence of 8 S from NUT charge shows that gravitomag- 
netic charge forces to shadow to get the shape of circle 
than ellipse. In the case of rotation, with the increase 
of black hole's angular momentum the shape of black 
hole shadow takes form of ellipse rather than circle. 
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4 PARTICLE CAPTURE CROSS SECTIONS 
FOR BLACK HOLE WITH NUT CHARGE 

In this section we will study the pure effect of NUT pa- 
rameter assuming that the the angular momentum of 
the black ho le is equal to zero. It has been shown in 
the paper of Abduiabbarov et al. ( 20081) that variables 
in the Hamilton- Jacobi equation for the particle motion 
around NUT black hole can be separated in the equa- 
torial plane. In the space-time metric ((lj we assume 
that the central object is non-rotating and particles are 
confined at the equator ial plane (a = 0, and 9 = n/2) 
It was first shown bv lZimmerman and Shahir (jl989t) 
for the spherical symmetric case (NUT spacetime) and 
later in the paper of ? for the axial symmetric case 
(Kerr-Taub-NUT spacetime) that the orbits of the test 
particles are confined to a cone with the opening angle 
9 given by cos 9 — 2£l/C. It also follows that in this 
case the equ ations of motion on the co ne depend on / 
only via I 2 ( j Abduiabbarov et ai1l2008l) . 

The main point is that the small value for the up- 
per limit for gravitomagnetic moment has been ob- 
tained by comparing theoretical results with experi- 
mental data a s (i) l/M < 10~ 24 from th e gravitational 
microlensing ( Rahvar and Habibil 12004) . (ii) l/M < 
1.5T0~ 18 fro m the interferometry experimen ts on ultra- 



cold atoms ( Morozova and Ahmedov 20091 ). (iii) and 
similar limit has been obtained fr om the experiments 



on Mach-Zehnder interferometer ( Kagrama nova et al 



2008) (here M is the total mass of central gravitating 
object). 

Due to the smallness of the gravitomagnetic charge 
let us consider the motion in the quasi-equatorial plane 
when the motion in 9 direction changes as 9 = it/2 + 
S9(t), where 88(t) is the term of first order in I, then it 
is easy to expand the trigonometric functions as sin 9 = 
1 - S9 2 (t)/2 + 0{59 4 (t)) and costf = 59{t) - 0(S9 3 (t)). 
Neglecting the small terms 0(S9 2 (t)), one can write 



2 4 6 8 10 

r/M 

Fig. 6 The radial dependence of the effective potential of 
radial motion of the massive particles for the different values 
of the gravitomagnetic charge: solid line for l/M = 0.1, 
dashed line for l/M = 0.5, and dot-dashed line for I = 0.9. 



motion of the massive particles for the different values 
of the gravitomagnetic charge is presented in Fig. [5] 

Assuming that the uncharged particle is moving 
slowly at infinity, i.e. £ ~ 1 one can easily rewrite 
the expression (TT9")) in the following form: 



R{p) = p 3 
where 

P = 



l-C)p 2 + C P 



r 



M 



IC 
~2 



C 



(21) 



Gravitational capture of the particle occurs for C < 
C CT . For the C = £ cr orbit spirals into a circular orbit, 
the radius of which is determined by the value of the 
multiple root of the polynomial (|21j) . i.e. discriminant 
of the later should vanish. Now it is easy to find the 
expression for £ cr as 



+ f, A = 



2Mr - I 2 , and X = and con- Ca = 16AP - QV 



sequently the equation of motion for the radial motion 
takes the following form 



^) =R(r)= [£ 2 -l-2U cS (r,l,C)]r i , (I!)) 



where £ and C is the energy and angular momentum of 
the particle per unit of its mass and the quantity 



16M 2 



(22) 



U eS (r,l,C) 



I 2 + Mr A£ 2 
E + 2E 2 " 



(20) 



can be interpreted as effective potential of the radial 
motion of the test particle at equatorial plane. The 
radial dependence of the effective potential of radial 



In the FigJ3 the dependence of C CT from dimension- 
less NUT parameter is presented. The dependence 
shows that the presence of the gravitomagnetic charge 
decreases the capture cross section for particles by black 
hole. 



5 Conclusion 

In this paper, we have studied the shadow of black 
hole with nonvanishing gravitomagnetic charge and an- 
alyzed how the shadow of the black hole will be dis- 
torted by the presence of the NUT parameter. From 
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l/M 



Fig. 7 The dependence of the critical angular momentum 
for capturing by central black hole from gravitomagnetic 
monopole momentum. 

the numerical results we have obtained that the NUT 
parameter forces to increase the size of the black hole 
shadow. The dependence of the distortion parameter S s 
from the NUT charge shows that the gravitomagnetic 
charge forces black hole's shadow to get the shape of cir- 
cle than ellipse. We have also studied the capture cross 
section for massive particles by black hole with non- 
vanishing gravitomagnetic charge and found its strong 
dependence from the NUT parameter. 
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Fig. 4 Silhouette of the shadow cast by a Kerr-Taub-NUT black hole situated at the origin of coordinates with inclination 
angle 9 — it/2, having a black hole angular momentum a and a NUT charge I. Upper row, left: a/M — 0.5, l/M — 0.1 
(solid line), l/M = 0.5 (dashed line), and l/M = 0.9 (dashed-dotted line). Upper row, right: a/M = 0.7, l/M = 0.1 (solid 
line), l/M = 0.5 (dashed line), and l/M = 0.9 (dashed-dotted line). Lower row, left: a/M = 0.8, l/M = 0.1 (solid line), 
l/M = 0.5 (dashed line), and l/M = 0.9 (dashed-dotted line). Lower row, right: a/M = 0.99, l/M = 0.1 (solid line), 
l/M — 0.5 (dashed line), and l/M = 0.9 (dashed-dotted line). The shadow corresponds to each curve and the region inside 
it. 
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Fig. 5 Observables R a and 8 a as functions of the NUT charge, corresponding to the shadow of a black hole situated at 
the origin of coordinates with inclination angle 9 — tt/2 and spin parameters a = 0.99. 
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